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ABSTRACT 

AU Microscopii’s debris disc is among the most famous and best-studied debris discs, and one of only two resolved debris discs 
around M stars. We perform in depth collisional modelling of the AU Mic disc including stellar radiative and corpuscular forces 
(stellar winds), aiming at a comprehensive understanding of the dust production and the dust and planetesimal dynamics in the 
system. Our models are compared to a suite of observational data for thermal and scattered light emission, ranging from the ALMA 
radial surface brightness profile at 1.3 mm to spatially resolved polarisation measurements in the visible. Most of the data are shown 
to be reproduced with dust production in a belt of planetesimals with an outer edge at around 40 au and subsequent inward transport 
of dust by stellar winds. A low dynamical excitation of the planetesimals with eccentricities up to 0.03 is preferred. The radial width 
of the planetesimal belt cannot be constrained tightly. Belts that are 5 au and 17 au wide, as well as a broad 44au-wide belt are 
consistent with observations. All models show surface density profiles increasing with distance from the star up to » 40 au, as inferred 
from observations. The best model is achieved by assuming a stellar mass loss rate that exceeds the solar one by a factor of 50. 
While the spectral energy distribution and the shape of the ALMA radial profile are well reproduced, the models deviate from the 
scattered light observations more strongly. The observations show a bluer disc colour and a lower degree of polarisation for projected 
distances < 40 au than predicted by the models. The problem may be mitigated by irregularly-shaped dust grains which have scattering 
properties different from the Mie spheres used in this work. From tests with a handful of selected dust materials, we derive a preference 
for mixtures of silicate, carbon, and ice of moderate porosity. We also address the origin of the unresolved central excess emission 
detected by ALMA and show that it cannot stem from an additional inner belt alone. Instead, it should derive, at least partly, from the 
chromosphere of the central star. 

Key words, circumstellar matter - stars: individual: AU Mic (GJ 803, HD 197481) - submillimetre: planetary systems - scattering - 
polarisation - methods: numerical 


1. Introduction 

Debris discs are remnants of the planet formation process and 
consist of (unobservable) planetesimals and collisionally replen¬ 
ished dust. They are ubiquitous around main-sequence stars 
with an incidence rate of about 20% f or FGK stars with ages 
from hundreds of Myr to several Gvr dEiroa et alJl2013l) . For 
10 - 800Myr-old A stars. ISu et al. (^ 06h auo t eda h igher de¬ 
tection rate of » 30%, whereas Thureau et al.l (|2014|) found a 
similar frequency as for FGK-type stars. The possible differ¬ 
ence in the debris disc occurrence rate between the FGK and 
A stars appears to be largely due to the observed age ranges 
and t he typical disc dec a y timescale for different spectral type s 
(e.g., iDecin et ^ l2003t IWvatt et ^ l2007t iKains et ^ l2()l ll) . 
The fre quency of debris discs around M stars remains contro¬ 
versia l ( Lestrade et al.ll20^ Gautier et al. 120071 : Forbrich et alJ 


for which several reasons could exist: (i) M dwarfs have low lu¬ 
minosities and are older on average than other stars so that their 
discs are faint and already mostly collisionally depleted, (ii) The 
blowout of dust around M stars may be favo ured since they 
likely po ssess strong stellar winds (e.g., recentlv IJohnstone et al.l 
l2015allbh . (iii) Owing to their low stellar mass, it is easy to strip 
planetesimals from an M-dwarf disc during close encounters 
with massive objects, (iv) Debris disc surveys, including stars 
of various spectral types, reveal generic temperatures for warm 
and cold disc components of ~ 190K and ~60K, respectiv ely 
jMorales et al.ll201 UlBallering et alj|2013l:l^en et al.ll2014l) . If 
this also holds for M stars, the M-type discs would be close to the 
star, meaning a faster dynamic timescale of the disc and a quick 
depletion of the dust reservoir. Around a dozen o f non-resolved 
M-star discs have been found in the recent ye a rs (iForbrich et aF 


I 2 OO 8 I: iPlavchan et al.ll2009r Lestrade et alJ 2009 . l2Q12h . D espite 

the h igh abundance of M dwarfs in the Galaxy (~ 80%, iLadal 
I 2 OO 6 I) . only a few discs around them have been detected so far. 


I 2 OO 8 : [^vchan et alj|2009t IChen et aT]|2014l) . [Theissen & West 


(|2014|) found a much larger number of old M dwarfs showing 
mid-IR excesses, interpretable as circumstellar dust emission. 
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However, due to the high fractional luminosities measured, the 
dust is more likely attributed to planetary impacts within the ter¬ 
restrial zone than to asteroid belt-like debris di scs. The discs 
aroun d the M stars AU Microscopii and GJ 581 (iLestrade et alJ 
12012|) are the only ones that are spatially reso lved thus far. 

_ Since the discovery of the AU Mic disc (iKalas et al.ll2004t 

lLiull2004l) it remains remarkable among the resolved debris sys¬ 
tems in many respects. The disc is seen edge-on with an impres¬ 
sive radi al extent of about 150 au and is resolved in the optical to 
nea^l R (iKrist et al.ll2005l: iMetchev et al.ll206^ iFitzgerald et al.l 
l2007h where it appears blue relative to the star. The surface 
brightness prohles from the optical/near-lR images show shal¬ 
low inner slopes at small projected distances but steepen sub¬ 
stantially beyond 35 au. In addition, asymmetries on small and 
large scales with several local brightness maxima and mi nima 
have been detected at stellocentric distances beyond 20 au (iLiul 
l2004tl&ist et al.ll2005HFitzgerald et aPUPOTh . In the latest high 
spatial r esolution STIS (Space T elescope Imaging Spectrograph) 
images, ISchneider et alJ (120141) found a distinct brightness en¬ 
hancement above the disc midplane on the south-east (SE) side 
at about 13au from the star. Furthermore, a disc warp is dis¬ 
cernible on the north-west (NW) side between 13 and 45 au, 
opposed to the SE brightness bump. There is also an NW-SE 
asymmetry, with the NW side brighter than the SE side in¬ 
terior to 20 au. These disc inhomogeneities hint at the exis¬ 
tence of planetary perturbers, causing radially localised struc¬ 
tures such as ri ngs, clumps, and gaps through planet-disc in - 
teractions (e.g.. lErtel et alJ l2012bt iNesvold & Kuchneiil2015h . 


However, no co nfirmation of planets in the AU Mic systern 
is found to da te (Neuhauser et al 2003t iMasciadri et alJ 120051: 
IMetchev et al.l l^05t [Hebb et al.l 2007h . Alternatively, clumpy 
disc structures may also be due to recent breakups of large plan- 
etesimals (e.g.. lKral et al.l[2()T^ . 

The disc has been also resolved at 1. 3 mm with SMA 
(Submillimetre Array), see IWilner et al.l (l2012h . and with 


ALMA (Atacama Larg e Millimetre/submillimetre Array), see 
iMacGregor et alJ (|2013|) . From the ALMA observations two dis¬ 


tinct emission components have been identihed: a dust belt that 
extends up to 40 au and a central emission peak that remains un¬ 
resolved. The dust belt shows an emission profile rising with the 
distance from the star, indicating a steep surface density slope. 
The central em ission peak is 6 times brighter than the stel¬ 
lar photosphere. IMacGregor et’aD (1201 3l) surmised this emission 
ste ms from an inner plane tesimal belt, located < 3 au from the 
star. ICranmer et al.l (1201 3h proposed it to come solely or partly 
from an active stellar corona. 

This great abundance of available observational data of the 
AU Mic disc provides valuable information from the smallest 
up to the largest dust particles, and further to planetesimals that 
produce the dust through collisional grinding. Constraints on the 
properties and spatial distribution of (sub)micron-sized particles 
can be inferred from the short-wavelength observations where 
the dust’s scattered light is detected. However in the AU Mic 
edge-on system, this is impeded by the degeneracy between the 
scattering properties of grains and their supposed spatial distri¬ 
bution, and requires a good knowledge of the phase and polari¬ 
sation function of small dust. The dynamics of small particles is 
strongly affected by radiative and corpuscular forces, resulting 
in radiative or stellar wind blowout and drag, depending on the 
grain and stellar properties. Observations at longer wavelengths 
show the thermal emission of larger grains and give insights into 
their spatial distribution. Especially, resolved images in the radio 
range reveal the locations of mm-sized, tracing the underlying 
parent bodies (planetesimals). 


Much work has been done in the past to chara cterise and 
model the AU Mic system. lAugereau & Beusll (l2006h performed 
a direct inversion of the visible and near-IR surface brightness 
prohles. From //-band observations a disc surface number den¬ 
sity was obtained that peaks around 35 au, close to the location 
of the break in the brightness prohle, which hints at a plan¬ 
etesimal belt around that distance. Furthermore, by htting the 
blue disc colour they constrained the grain size distribution and 
found that grains smaller th an 1 pm are mandatory to e xplain the 
observations. The study of IStrubbe & ChialigI (l2006h included 
the dynamics of grains for the hrst time. There, a narrow birth 
ring of planetesimals at about 40 au was assumed. Micron-sized 
dust grains are produced through mutual collisions of the birth¬ 
ring objects. Grains with sizes < 1 pm that are smaller than the 
radiative and corpuscular blowout limit, are expelled from the 
system by direct radiation and stellar wind pressure. Slightly 
larger grains are either barely bound and launched into eccen¬ 
tric orbits, hence forming a halo, or are transported into the 
inner regions of the system by Poynting-Robertson and stellar 
wind drag. The largest grains foll ow the dynamics o f their par¬ 
ents and remain in the birth ring. IStrubbe & Chi an a concluded 
that the disc is collision dominated with an inner part devoid 
of submi cron-sized grains (co nsistent with polarisation measure¬ 
ments of iGraham et al.ll2007l) . and small grains mainly populate 
the outer p art of the disc, causing the blue colour of the scat¬ 
tered light. IFitzgerald et al.l (120071) attempted to simultaneously 
reproduce the scattered light prohles, the degree of polarisation, 
and the spectral energy distribution (SED) with a two-zone disc, 
using power-law descriptions for the radial and size disPibu- 
tion of the dust in each zone. They showed that the hrst zone 
(35 - 40 au) of large particles, representative for a planetesimal 
ring, mainly accounts for the long-wavelength thermal emission, 
whereas the second zone (40 - 300 au), composed of small parti¬ 
cles, well reproduces the scattered light measurements. Recently, 
a dust halo was also detected in far-IR resolved Herschel and 
JCMT (James Clerk Maxwell Telescope) images by Matthews et 
al. (in prep.). Their best ht s upports the narrow birth-ring model 
of IStrubbe & Chian^ (l2006l) . but is also consistent with an ex¬ 
tended pla netesimal ring from 8 - 4 0 au, according to the mo del 
explored in IWilner et all (1201 2l) and iMacGregor et aP (l2013l) . 


In this study, we undertake collisional modelling to hnd a 
generic size and radial distribution of the dust in the AU Mic 
disc. We aim at a comprehensive understanding of the dust pro¬ 
duction and the dust and planetesimal dynamics in this sys¬ 
tem. Eor the hrst time, we combine constraints from scattered 
light and mm wavelength observations and search for a self- 
consistent collisional model that explains all these data. We as¬ 
sume an axisymmetric disc, and therefore do not account for 
the formation of substructures and asymmetries observed in the 
AU Mic disc. Section |2] gives an overview of the observational 
data and presents our re-reduction of the ALMA data, resulting 
in a star-Hdisc Hux density and a radial surface brightness prohle 
at 1.3mm. Section [3] explains our collisional model. Section]?] 
shows the modelling of the extended, ss40au-wide emission 
zone, seen in the ALMA image, which we refer to as resolved 
outer disc. Section ]5] addresses a possible origin of the unre¬ 
solved central emission. In Sect.]6] we compare our results from 
collisional modelling with a multidimensional power-law htting 
approach. Conclusions are drawn in Sect. ]7] 
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2. Data used 

2. 1. Stellar properties 

AU Mic (GJ 8 03, HD 197481) is an Ml V dwt^f at a dis¬ 
tance of 9.9 pc (iPerrvman et alj[l997l : Ivan Leeuwenll2007h . The 
star is a brigh t X-ray and UV emitt er and shows strong flaring 
activity ( e.g.. iRobinson et al.l 20W). We used the photosphere 
model of lAugereau & Beusj ( 2006 ) with an effective tempera¬ 
ture of 3700 K, a luminosity of 0.09 Lq, and a surface gravity 
of log(g) = 4.5 (CGS). We assumed a stellar mass of 0.5 Mq, 
motivated by the me an of the wide mass range given in the lit¬ 
eratur e (0.3 - 0.6 MfT^. lPlavchan et alJ2009t[Houdebine & Dovlel 
Il994l) . Note that this assumption is roughly consistent with the 
photosphere model used that predicts » 0.6 Mq. 

AU Mic belongs to the p Pic moving group (BPMG) and its 
age coincides with the BPMG age by assuming a coeval origin of 
all group membe rs. Through the identific atio n of the lithium de- 
pletion boundarv. lBinks & JeffriesI (1201 4t) and iMalo et al.l (l2014t) 
found BPMG ages of 21 + 4Myr and 26 + 3 Myr, r espec¬ 
tive ly. This agrees w i th tra ceback ages from iMakarovI (l2007h 
and iMamaiek & Bell (l2014h . The latter used revised Hipparcos 
astrometry data and also derived an isochronal age of 22 + 3 Myr. 


2.2. ALMA observations 


ALMA observations of AU Mic in Band 6 (1.3 mm) have been 
carried out as part of the cycle 0 early science observations 
in the context of the projects 2011.0.00142.S (PI; D. Wilner) 
and 2011.0.00274.S (PI: S. Ertel). The data considered in the 
present work were taken on 16 June 2012 and are in agree¬ 
ment with the iMacGregor et al.l ( l2013h ’s SB-4 observations. A 
total of 20 operational 12-m antennae were used, spanning base¬ 
lines of 21 - 402 m with an effective angular resolution (FWHM 
of the reconstructed dirty beam) of 0.69" x 0.79" and an ef¬ 
fective field of view of A/D » 22" (with T the observing 
wavele ngth and D the single dish diameter). IMacGregor et alJ 
(1201 3l) presented a detailed description of the observations 
and data reduction. Our re-reduction of the data was carried 
out in the Common A stronomy Software Applications (CASA, 
McMullin et aO l2007h . We followed the approach used by 


MacGregor et alJ(l2013h and came to consistent conclusions. 

The edge-on disc has a radial exte nt of about 4" (sa 40au ) 
and a position angle of 128.41 + 0.13° (IMacGregor et alJl2013h . 
It is unresolved in its vertical direction. We refrain from pre¬ 
senting the ALMA image as it would not provide inf ormation 
beyond that shown in Fig. 1 of iMacGregor et akl (l2013l) . In order 
to prepare the data for our modelling, we extracted a radial pro¬ 
file along the disc major axis. To this end, we largely adopted 
the approach used for Herschel data obtained in the context 
of the Op en Time Key Programme DUst around NEarby Star s 
(DUNES. lLdhne et al.l20ialEiroa et alJ2013l:lErtel et alJ20l4 . 
We first converted the flux units in the image from Jy/beam to 
Jy/pixel by multiplying them with a factor of 0.0185 derived by 
integrating the core of the reconstructed dirty beam. Then, we 
rebinned the image with an original pixel scale of 0.1 "/pixel to 
a ten times smaller pixel scale using a cubic spline interpola¬ 
tion and rotated it by an angle of 38.41° clockwise in order to 
align the disc major axis with the image x-axis. From this im¬ 
age, we measured the radial profile along the disc major axis by 
integrating over 11x11 subpixel-wide boxes centred on the disc 
midplane and spaced by 74 subpixels (0.74", about the FWHM 
of the reconstructed dirty beam, i.e. the resolution element of 
the image) left and right of the star, assumed to be located at the 


brightest subpixel. We averaged the SE and NW sides of the pro¬ 
file to fit an axisymmetric model to the data. Uncertainties were 
derived as a quadratic sum of the difference between the two 
sides of the profile and the background fluctuation measured in 
regions of the image where no flux is expected but where the 
sensitivity is comparable to the image centre. As we were only 
interested in the dust distribution in the outer disc, we ignored 
the profile point at the image centre, which is affected by the 
inner, unresolved component. By extrapolating the disc profile 
from the outer points to the disc centre and assuming that the 
flux measured above the extrapolated value stems from the inner 
component (star and additional emission), we estimated a flux 
density of this c omponent of 0.29 + 0.05 mJy, consistent with 
0.36+0.07 mJy in fMacGregor et ai](l2013l) . We estimated the un¬ 
certainty of this measurement by adding in quadrature a typical 
uncertainty in the outer points (~ 15% ) and an absolute photo - 
metric calibration uncertainty of 10% (IMacGregor et al.ll2013l) . 
The stellar contribution to this flux de nsity is 0.04 mJy, obtained 
from the stellar photosphere spectrum (lAugereau & Beus'3l2006h 
extrapolated to the ALMA wavelength using the Rayleigh-Jeans 
law. However, this might significantly underestimate the stellar 
flux and the whole central compon ent might originate from stel¬ 
lar emission (ICranmer et alJl201^ . 

Deviating from th e complex fitting approach used by 
IMacGregor et al.l (1201 3h to derive the disc parameters (including 
the total flux density), we performed photometry of the disc by 
integrating the flux in the rotated image in a box of 101 x 21 
native image pixels. The uncertainty was estimated from the 
scatter of the flux measured on eight positions above and be¬ 
low the disc, where no emission is expected. An additional cal¬ 
ibratio n uncertainty of 10%, consistent with IMacGregor et al.l 
(1201 3h . was added in quadrature. We found a total flux density 
of 8. 75 ± 0.98 mJy, w h ich is consistent with 7.46 ± 0.76 mJy 
from IMacGregor et akl (1201 3l) . by summing both disc compo¬ 
nents and including calibration uncertainty. For the modelling, 
we subtracted the flux of the inner component estimated above, 
so that we are left only with the flux of the disc and the star as 
estimated using our photosphere spectrum. 

Note that the exact flux of the inner component has no signif¬ 
icant impact on our profile or flux measurement, since we ignore 
the inner profile point affected by this component and because 
its contribution to the total flux density is only 0.2cr. 


2.3. Auxiliary data 


In addition to the ALMA data, we considered for modelling 
a variety of photometric data fro m earlier work . Furth ermore, 

1201 Oh PACS 


we 


included archival Herschel (iPilbratt et al 


(Phot odetector Array Camera and Spectrometer, 


Poglitsch et al.l 


201 oh and SPIRE ( Spectral and Photome tric Imaging Receiver, 
Griffin et'aP 1201 Ot ISwinvard et al.1 1201 Oh data, obtained in the 
context of the debris disc Guaranteed Time Observations pro¬ 
gramme (PI: G. Olofsson, see Matthews et al., in prep.). These 
data were trea ted following the i/eric/zel/DUNES data reduc¬ 
tion strategy (lEiroa et al. l2013h. HIPE {Herschel Interactive 
Processing Environment, ' lOttll2010h version 11 and the calibra¬ 
tion plan versions v56 and vll for PACS and SPIRE, respec¬ 
tively, were used. PACS flux densities at 70 p.m and 160 pm 
were measured with an aperture of 20" and annuli for the noise 
measurement of 30 - 40" and 40 - 60", respectively. The SPIRE 
flux density at 250 pm was measured with an aperture of 30" 
and the noise was estimated in a region of the image where no 
emission from the disc is expected. A nearby source is visible 
in all Herschel images, peaking at 250 pm. It is well separated 
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from the disc at wavelengths up to 250 iim and easily removed 
by a point source subtraction. The two sources are not well sep¬ 
arated at longer wavelengths. SPIRE flux de nsities at 350 lam 
and 5 00 p.m were measured using SExtractor (iBertin & ArnoutsI 
[l996h and a point source approximation (PSE photometry), since 
the source is unresolved at these wavelengths. Aperture correc¬ 
tions were applied to the measured fluxes. Instrument calibration 
uncertainties of 7% for PACS and 6% for SPIRE were added in 
quadrature to the measured uncertainty. Note that a more de¬ 
tailed reduction of the Herschel data will be given in Matthews 
et al. (in prep.). Minor differences to the data used here exist, but 
all within 2cr (B. Matthews, priv. comm.). The differences have 
no impact on the results presented in this work. 

We did not consider the spatially resolved information from 
the Herschel images since the resolution is much lower than that 
of our ALMA image (e.g., FWHM — 5.8" for HerschelfPACS 
at 70 pm vs. FWHM - 0.74" for ALMA Band 6). However, 
we inspected the PACS data for signs of bright, warm central 
emission in the form of a s ignificantly incre asing extent of the 
disc from 70 pm to 160 pm (lErtel et al.ll2014lf . Such a behaviour 
is not visible suggesting that the inner component seen in the 
ALMA images has indeed no o r only a minor conPi bution to 
the fluxes at shorter wavelengths. iDoering et al.l (l2005l) reported 
a detection of the innermost disc regions up to a; 16au in di¬ 
ameter via A-band imaging at Gemini South. However, only a 
low-level surface brightness extension along the direction of the 
edge-on disc was found that has not been taken into account in 
our study. Eurthermore, we did not consider Spifzer/MIPS im¬ 
ages due to the low angular resolution (e.g., the resolution of 
Spitzer at 24 pm is nearly as low as that of Herschel at 70 pm). 
Any confusion, e.g. from a background galaxy, would here be 
difficult to disentangle from the source. 

Assuming the flux density to be proportional to at long 
wavelengths, the mean specttal index beyond 250 pm amounts 
to (5 = 1.7, close to the Rayleigh-Jeans (RJ) slope of a black- 
body radiator (cJrj = 2). We fitted a blackbody curve to the far- 
IR excess in order to derive spectral slopes at all wavelengths. 
A temperature of 50 K fits the data very we ll, consistent with 
iRebull et alJ(l2008t) and IPlavchan et al.](l2009l) . Interestingly, this 
is close to the generic temperature found for the o uter cold com¬ 
ponents in many two-component disc systems jMorales_e^^ 
20n; iBallering et al.l 1201 3t IChen et ^ l2014t iPawellek et al. 
20141) . The temperature that was found for the inner disc, 190 K, 
corresponds to a very small distance of 0.7 au for a putative 
warm component in the AU Mic system. This is well below the 
resolution limit and leaves room for further speculations about 
the presence of an inner unresolved disc component. We took 
the nearest tabulated values for the colour corrections listed in 
the literature for the instruments used and the specPal slopes 
of the source (star and disc) derived from this fit. The resulting 
colour corrected flux densities and references for both the mea¬ 
surement and the colour correction are listed in Table [T] 

2.4. Scattered light 

The AU Mic disc is 60% brighter i n B-band than in /-band rel¬ 
ative to the star (iKrist et al.ll2005h . A potential explanation is 
a significant amount of small grains th a t scatt er more light at 
shorter wavelengths. lAugereau & Beus3 (l2006l) constrained the 
grain size disPibution by fitting the blue disc colour. They con¬ 
cluded that grains down to ~ 0.1 pm are necessary to explain 
the observati ons when astronomic al silicate or graphite grains 
are assumed. iGraham et akl (l2007h obtained polarisation maps 
in U-band using the Hubble Space Telescope (HST) Advanced 


Table 1. Colour-corrected photometry of AU Mic. 


Wavelength 

[pm] 

Flux density 
[mJy] 

Telescope/ 

Instrument 

Ref. 

Elux 

Ref. 

Col. corr. 

3.35 

4820±119 

WISE 

C12, 

WIO 

WIO 

4.6 

4260±119 

WISE 

C12, 

WIO 

WIO 

5.8 

2200 ± 55 

SpitzerflRAC 

C05 

IH 

8 

1270 ± 20 

SpitzerflRAC 

C05 

IH 

9 

912 ±22 

Akari 

no 

IH 

11.6 

543 ± 60 

WISE 

C12 

WIO 

12 

517 ±30 

IRAS 

M90 

B88 

18 

246 ± 36 

Akari 

no 

IH 

22.1 

183 ±21 

WISE 

C12 

WIO 

24 

158 ±3 

Spitzer/MlPS 

P09 

IH 

25 

244 ± 63 

IRAS 

M90 

B88 

60 

269 ± 46 

IRAS 

M90 

B88 

70 

227 ± 27 

SpilzerfMlPS 

P09 

IH 

70 

231 ± 16 

HerschelfPACS 

HSA 

IH 

100 

680 ±149 

IRAS 

M90 

B88 

160 

172 ±21 

SpitzerfMlPS 

R08 

IH 

160 

243 ± 17 

HerschelfPACS 

HSA 

IH 

250 

134 ±8 

HerscheljSPtKE 

HSA 

IH 

350 

72 ±21 

CSO/SHARCII 

C05 


350 

84.4 ±5.4 

HerschellSPWE. 

HSA 

IH 

450 

85 ±42 

JCMT/SCUBA 

L04 


500 

47.6 ±3.8 

HerscheljSPlRE 

HSA 

IH 

850 

14.4 ± 1.8 

JCMT/SCUBA 

L04 


1300 

8.5 ±2.0 

SMA 

W12 


1300 

8.75 ±0.98 

ALMA 

M13 



B88: iReichman etZI ( fIMsh . C05: IChen et all ( l200l . 

iDE 


Refe rences._ _ _ _ _ _ 

C12: ICutri & et all ll2012h . HSA: Herschel Science Archive (Matthews 
et al., in prep.). IH: corresponding instrument hand book. 110: 
llshihara et al] OOl^. L0 4 :[Liu et all d^Ol), M90:lMoshir et al.l 

' jgf). 


Ml3: iMacGregor et all (l2013h. PQ9: iPlavchan et all (l2009l). R08: 


IRebull et all j2008h . WIO: IWright et all ( |201(]|) . '^2: IWilneretan 
(l2012h . The ALMA and Herschel data were re-analysed in the present 
work. 


Camera for Surveys. They measured the degree of linear po¬ 
larisation to be steeply rising from 5 to 40% within a distance 
of 80 au from the star. The light is polarised perpendicular to 
the disc plane that agrees with the expected scattering behaviour 
by small grains with size parameters x < I (equivalent to typ¬ 
ical grain sizes < 1 pm). To model th e V-band scattered l ight 
intensity and the degree of polarisation. lGraham et all (l2007h as¬ 
sumed uniformly disPibuted dust between two distances from 
the star and adopted the empirical iHenvev & GreensteinI (Il94lh 
phase function in combination with a parameterized polarisation 
function. In their best-fit model, the particles exhibit sPong for¬ 
ward scattering and the innerm o st 40 au of the disc are devoid 
of gr ains. Later on , JShen et al.l (l2009h showed that the results 
from iGraham et al.l (120071) can be reproduced by a distribution 
of sphere cluster aggregates with size disPibution index -3.5, 
volume-equivalent radii between 0.1 and 0.4 pm, and a porosity 
of 60%. 


In this work, our models are compa red with the opti c al and 
nea r-IR surface brightness profiles from lEitzgerald et^ (l2007h 
and lSchneider et ah ( 20141) an d with t he measured degree of po¬ 
larisation from Graham et al.l (l2007h . As our modelling deals 
with axisymmetric discs, all observed profiles had been averaged 
over the SE and NW disc sides in order to remove the observed 
surface brightness anomalies. 
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2.5. Gas in the disc 


Some studies searched for circumstellar gas in the AU Mic disc. 
From a non-detection of submillimetre CO emission Liu et alJ 
(l2004ll inferred an H 2 mass of < 1.3 Mm. iRoberge et ^ ( 2005h 
reduced the H 2 mass l imit to 0.07 Mm due to a non-detection of 
far-UV H 2 absorption. iFrance et alJ (l2007h analysed fluorescent 
H 2 emission and found a total gas mass between 4 x 10 
and 6 x 10“^ Me, but stated that this detection might also come 
from a c loud that extends beyond th e disc. From X-ray spec¬ 
troscopy ISchneider & Schmiti (1201 Oh reported a maximum H 
column density of ~ 10'® cm“^, which is about five times higher 
than the interstellar value. However, since observational data 
point towards a very low gas content in the AU Mic disc, we 
neglect the influence of any gas in our investigations. 


3. Collisional modelling 


Collisional modelling with our code ACE (Analysis of 
Collisional Evolution) and its application to debris disc systems 


is described extensively in previous papers ( 

Krivov et al. 120061 

2008 

Muller et al. 

2010 

Reidemeister et al. 

2011; iLohne et al.l 

2012 

Krivov et al. 

2013 

Schiippler et al.ll20l4 among others). 


The code assumes a disc of planetesimals orbiting a central star, 
simulates their collisional grinding, and follows the collisional 
and dynamical evolution of the system under a variety of physi¬ 
cal effects. ACE works with a three-dimensional logarithmic grid 
for objects mass, pericentric distance, and eccentricity. The sim¬ 
ulations consider a number of parameters, including those that 
define the initial distribution of planetesimals, the strength of 
particle transport, and material properties such as the critical 
fragmentation energy. The following sections detail the parame¬ 
ters chosen in our modelling. 


mass loss of M* = 300 Mq. This high wind strength means that 
the direct stellar wind pressure is larger than the radiation pres¬ 
sure which causes an efficient transport of small dust into outer 
disc regions. 

In our modelling, we considered moderate and strong wind 
strengths by assuming SOM© and 300M©- Throughout this pa- 
per, we refer to bo t h case s as SW50 and SW300. According to 
lAugereau & Beus3 (l2006h . SW50 is a typical model where the 
flares are ignored, while SW300 is an attempt to account for the 
impact of the episodic stellar activity. In all simulations involv¬ 
ing stellar winds, we used Eq. (7) of iGustafsonI d 19941) for the 
total stellar wind force a cting on a particle and a ssumed a wind 
speed Vsw = 400 km s ' dStrubbe & Chiang||2006l) . 


3.2. Grain materiai 


iGraham et al.l d2007h and lEitzgerald et al.l d2007h found evi¬ 
dence for porous aggregates containing silicate, carbon, and 
ice. Guided by these res ults, we used porous grains composed 
of as tronomical silicate dDrainell2003h . ACH2 amorphous car - 
bon dZubko et al.lll996l) . and water ice (Lj_&_Greenbergl [l998l) . 
Since the ACH2 refractive data from IZubko et alJ ( 19961) termi¬ 
nate at 4 « 1 mm and do not cover the full wavelength range 
for all available photometric measurements, we extrapolated 
the real and imaginary parts up to 4 = 2 mm by the function 
Cl (4/1 mm)'^^. The extrapolation gives (ci,C 2 ) = (16.2,0.27) 
and (ci, C 2 ) = (5.6, -0.26) for real and imaginary parts, respec¬ 
tively. 

We generated different porous silicate-carbon-ice mixtures 
and calcul ated the refr a ctive indexes of the composites by ap¬ 
plying the iBruggenianl d 19351) mixing rule. The materials con¬ 
sidered in this study are given in Table |2] Eor brevity, we use 
the identifiers Ml to M5 when referring to the different compo¬ 
sitions. 


3.1. Steiiar winds and mass ioss rate 


Stellar winds (SW) are expected to play a crucial role in discs 
around late-type stars such as AU Mic where they are typ¬ 
ically stronger than for early-type s tars (e.g., Plavchanet_a]J 
2005L 200^ Strubbe & ChiangllioO^ iReiderneister et al.ll201 It 
Schiippler et al. 2014i) . The total stellar wind force that acts on 


the circumstellar dust can be decomposed to direct stellar wind 
pressure, which is directed radially outwards, and stellar wind 
drag, which works like a headwind on the particles that causes 
them to gradually lose orbital e nergy and angular momentum, 
and therefore, to spiral inwards (iBurns et alJll979l) . Depending 
on the stellar wind strengths and the size of the particles they 
are interacting with, the result could be either a grain transport 
towards the star or a blowout of smaller grains from the system. 

Previous studies proposed various stellar winds strengths 
for the AU Mic system, expressed in terms of multiples 
of the solar mas s loss rate, Mq = 2 x IO^'^'Mq yr“'. 

IStrubbe & Chianj (l2006l) constrained the stellar mass loss rate 
to be M* < IOMq which was in good agreement with the ob¬ 
served SED and the V - H colour profile of the disc. This 
relatively weak wind strength means that the disc is rather 
dominated by destructive grain- grain collisions than transport 
processes. IPlavchan et alJ (120091) estimated <50 Mq from 
consid erations of Spitzer/IRS observations. lAugereau & Beus3 
(I 2 OO 6 I) took into account the active nature of AU Mic, and found 
M* 50 Mq during quiescent states and M* 2500 Mq during 
flares. The star spends about 10% of the tim e in flare phases wit h 
typical flare durations of several minutes (iPagano et al.ll2000l) . 
Averaging over quiescent and flare phases yields a mean stellar 


Table 2. Materials used for modelling. 


Name 

Composition 

Density 
[g cm-3] 

Blowout size [pm] 

SW50 SW300 

Ml 

sil33-l-car33-(-vac33 

1.78 

0.07 

0.48 

M2 

sil504-car50 

2.68 

0.04 

0.35 

M3 

sil25 -l-car25-l-ice25 vac25 

1.64 

0.07 

0.51 

M4 

sil 10-t-car 10-l-ice40d-vac40 

1.02 

0.11 

0.71 

M5 

sill5-l-carl5-l-vac70 

0.80 

0.16 

0.93 


Notes. Numbers after material names denote their volume fraction in 
percent, e.g. sil25-l-car25+ice25+vac25 = 25% silicates + 25% carbon 
-I- 25% ice -t- 25% vacuum. The bulk densities of the composites are 
given by = 2p/ tr,, where cr, denotes the volume fraction of the 
species, and p, their densities (psn = 3.5gcm“^, pcar = 1.85gcm“^, 
Picc = 1.2gcm“^). Blowout sizes are grain radii for which j8 equals 0.5 
(see text). 


We assumed spherical dust grains. Each grain has a 
certain radius s which is denoted as its size. Optical 
properties of the grains were determined by means of 
Mie theory dBohren & Huffman! 119831) . Eigure [T] shows 
P = A-ad + Aw = Frad/Fgrav + Fsw/Fgrav, the ratio of direct radi- 
ation and direct st ellar wind pressure forces to stellar gravity 
(iBurns et al.lll979t) as a function of grain size. A detailed de¬ 
scription of Aw includes a dependency on the distance from 
the star, r, due to radial variations of the stellar wind velocity. 
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Ch. Schiippler et al.: Collisional modelling of the AU Microscopii debris disc 



grain size, s [pm] 


where (2 d,s, bs) and (Qo.g, bg) are two pairs of constants, char¬ 
acterising the strength and gravity regime, respectively, and Vimp 
denotes the impact speed. is the impact energy per unit tar¬ 
get mass necessary to disperse the target into fragments where 
the largest fragment contains at most half of the original tar¬ 
get mass. For objects smaller than 100 m, 2* is mainly deter¬ 
mined by the material strength, while for larger sizes the grav¬ 
itational binding energy dominates. The material strength also 
depends on the material composition, but this dependence is 
still insufficiently constrained by laboratory work. Therefore, 
we used the same 2* for all materials given in Table and 
assumed 2 n s = O.n.v - ^ x lO ^ergg^*, - -0.37, bg = 1.38, 
and K - 0.5 (iLdhne et al.ll2012i and references therein). 


3.3. Setup 


Fig. 1. yS-ratio for different assumptions of AU Mic’s stellar mass 
loss rate M*. The material consists of 33% carbon, 33% silicate, 
and 33% vacuum (Ml). The horizontal dashed line depicts /? = 
0.5. Particles with fi > 0.5 are ejected from the system owing to 
radiation and stellar wind pressure. 


Howe ver, this effect is weak for r > 15 au (lAugereau & Beusl 
|2006|) and we considered jSs* independent of r. Further, we as¬ 
sumed the reflection of wind particles on the surface of the 
dust gra ins that means a n efficiency of momentum coupling of 
Cd = 2 (lGustafsonlll994 . 

The orbits of the collisionally-produced dust grains become 
more and more eccentric with increasing [i. For fi > 0.5, the 
grains get eccentricities > 1 and are expelled from the system 
on hyperbolic orbits if they are released from bodies on circu¬ 
lar orbits. As a measure for the blowout size, we define ibiow = 
s{j3 - 0.5). However, when particles are released from eccentric 
orbits, the blowout occurs bet ween^ = (l-e)/2 at the periastron 
and (1 -H e)/2 at the apastron (iBurns et al.lll979l) . where e is the 
typical eccentricity of the parent bodies that eject the dust. ACE 
computes the orbital elements of the fragments depending on 
the masses, semi-major axes, eccentricities, and j6’s of the col¬ 
liding bodies, and thus, automatically accounts f or bound and 
unbou nd grains launched from non-circular orbits dKrivov et alJ 
1 20061) . Therefore, the minimum size of bound grains in the de¬ 
bris disc system is slightly smaller than Sbiow Nevertheless, we 
found this effect to be small and Sbiow can be considered a repre¬ 
sentative value for the minimum grain size. 

In the case of weak or no stellar winds, fi is smaller than the 
blowout limit 0.5 for all grain sizes (see Fig. [1] M* = 0). With 
increasing stellar wind strength the dust particles experience a 
stronger radial wi nd pressure. Assum ing a radiation pressure ef¬ 
ficiency 2pr = 1 (iBurns et alJll979h . yields Fsw/Fiad ~ 0.5 for 
SW50 and Fsw/Frad ~ 3 for SW300. As a result, if the stellar 
winds are sufficiently strong, the jS-ratio that would otherwise 
be smaller than 0.5 can exceed this critical value at the small¬ 
est dust sizes. This defines a blowout size, terminating the lower 
end of the size distribution (Table |2]i. Note that Fig. [1] shows fi 
for Ml exemplarily. The other materials given in Table |2] follow 
a similar dependence of fi on size. 

The mechanical str ength of the mater i als wa s defined by a 
modified version of the iBenz & Asphaud d 19991) description of 
the critical energy for fragmentation and dispersal. 


Q 


★ _ 

D ” 



+ Go.g 



^imp 

3 km s-‘ 


( 1 ) 


Table [3]lists all ACE simulations discussed in this paper. We dis¬ 
tinguish between three disc types, named by the width and loca¬ 
tion of their planetesimal belt (PB): (i) a narrow outer PB; (ii) a 
radially extended PB, and (iii) a narrow inner PB. Configurations 
(i) and (ii) describe models for the outer resolved disc, whereas 
(iii) corresponds to a model for the central emission seen by 
ALMA. In all simulations, we assumed the discs to be com¬ 
posed of planetesimals up to at least 10 m in radius with an ini¬ 
tial size distribution oc . The planetesimals are large enough 
to ensure that their collisional lifetime is longer than the time 
neede d to reach a quasi-s teady state material distribution at dust 
sizes (iLohne et alJ 120081) . Thus, the results of the simulations 
are independent of the maximum size of the objects. The code 
also accounts for cratering and bouncing collisions where the 
bulk of one or both colliders stays intact but fragments are re¬ 
leased. For each collision, the distribution of fragments is as¬ 
sumed to follow a power law with a mass distribution index 
-1.83 down to the smallest objects represented in the grid. The 
planetesimal discs were assumed to be stirred to a maximum 
eccentricity, emax- The bodies in the disc were assumed to be 
uniformly distributed within a semi-opening angle /max = eniax/2 
(energy equipartition). After reaching a quasi-steady state, each 
simulation was evolved until the dust mass had reduced suffi¬ 
ciently to reproduce the observed excess emission. In general, 
the simulation time, fsim, does not coincide with the physical age, 
fphys ~ 20 - 30 Myr, of the AU Mic system. For all models ex¬ 
cept the one with planetesimals between 1 and 45 au, fsim was 2 
to 4 times shorter than fphys. These times are more than enough to 
bring the model discs to collisional equilibrium. The initial con¬ 
ditions of our model discs are, however, fiducial. Different com¬ 
binations of total initial mass and initial size distribution can lead 
to the same observed dust configuration, though after different 
times. An estimate for the minimal initial disc mass necessary to 
provide the observed amount of dust after fphys can be obtained 
by assuming that the collisional lifetimes, Loi, of the largest bod¬ 
ies equal the system age, enabling them to sustain sufficient dust 
production over that time. The required sizes for the largest bod¬ 
ies, typically exceed those in our simulation runs. We ex¬ 
trapolated the size distributions and collisional lifetimes to de¬ 
rive such that fcoK^max) ~ ^physj then, determined the 
minimum initial disc masses Mmin for particle sizes up to 
The power-law distributions in the m-to-km size range facilitate 
that extrapolation. The results are listed in Table[3] These masses 
reflect the minimum disc masses necessary to achieve simulation 
times that would approach fphys. Larger planetesimals and poten¬ 
tial exo-Plutos cannot be assessed this way because their direct 
contributions to the collisional cascade are negligible. 
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Table 3. Description of ACE runs. 


Disc type 

Extent" 

r-distr.'’' 


*^max 


Material" 

Mf 

^max 

M . " 

mill 

Grid resol." 

Sect.’“' 


[au] 


[Mg] 


[rad] 


[Me] 

[m] 

[Me] 






0 




1.5 X 10“* 

20 

0.02 




37.5-42.5 

oc r® 

50 

0.03 

0.015 

Ml 

1.7 X 10“' 

100 

0.05 

ml,el,pl 

I4d1 




300 




2.1 X 10“' 

250 

0.1 









M2 

2.4 X 10“' 

100 

0.06 



narrow outer PB 

37.5-42.5 

oc r® 

50 

0.03 

0.015 

M3 

1.5 X 10“' 

100 

0.05 

ml,el,pl 








M4 

1.4 X 10“' 

200 

0.05 




37.5-42.5 

oc r® 

50 

0.01 

0.005 

Ml 

1.7 X 10“' 

30 

0.05 

ml,el,pl 

1431 





0.1 

0.05 


6.4 X 10““* 

100 

0.03 


25.5-42.5 

oc r® 

50 

0.03 

0.015 


1.7 X 10“' 

100 

0.08 

ml,el,pl 

EH 

extended outer PB 

OC 

Ml 


1.0-45 





2.1 X 10“' 

100 

0.2 

m2,el,pl 







Ml 

6.2 X 10“® 

1 X 10' 

0.02 



narrow inner PB 

2.8-3.2 

ocr® 

50 

0.03 

0.015 

M3 

5.4 X 10“® 

4x 10' 

0.02 

ml,el,p2 

m 







M5 

4.5 X 10“® 

4x 10' 

0.02 




Notes. “'Radial extent of the planetesimal belt; '’'Specification of the radial distribution of planetesimals; “'Stellar mass loss rate in units of the solar 
value; “''Maximum eccentricity of the planetesimal orbits; “'Semi-opening angle of the planetesimal disc; ''Name of dust material used (symbol 
explanation in TableO; ^'Derived dust mass (grains up to 1 mm in radius); '’'Object size for which the collisional lifetime equals the physical age 
of the AU Mic system; “'Estimated minimum total disc mass (grains up to size i’Jnjy,); "Specification of grid: 37 bins in mass from 5 x 10“'* to 
1 X lO'® g (ml), 40 bins in mass from 8 x 10“'* to 1 x lO" g (m2), 22 bins in eccentricity from 0.001 to 2 (e2), 30 bins in pericentric distance from 
1 to 60 au (pi), 20 bins in pericentric distance from 0.1 to 7 au (p2); '“'Section number in this work where the runs are addressed. 


4. Resolved outer disc 

First, we seek to explain the available observational data on the 
resolved outer AU Mic disc, including the SED, the ALMA ra¬ 
dial brightness profile, and scattered light measurements. In the 
following, we present several modelling attempts, illustrating the 
influence of different parameters and their compatibility with the 
data. We start with an investigation of the stellar wind strength 
to And an appropriate reference run (Sect. |4~T]) . In a subsequent 
analysis (Sects. l4~2l 14.31 and l4.41 i. we vary other initial disc pa¬ 
rameters of this reference run in order to more tightly constrain 
the parameters and reveal degeneracies among them. 


4. 1. Stellar wind strength 

In a first step, we sought after a reasonable assumption for 
the stellar wind strength and carried out the runs SWO, SW50, 
and SW300. SWO does not assume any stellar wind activity at 
all (M* = 0). Here, the transport of particles is controlled by 
Poynting-Robertson drag only. In SW50 and SW300, moderate 
and extreme stellar winds are considered, given by 50 times and 
300 times the solar strength, respectively (see Sect. [3T|for justi¬ 
fication of these values). 

For other disc parameters, we took reasonable values that 
have been prefer red by previous studies. I nspired by the birth¬ 
ring scenario of IStrubbe & Chianel (l2006l) . we assumed a ra¬ 
dially narrow planetesimal belt centred at a stellocentric dis¬ 
tance rpB = 40 au, having the full width ArpB = 5 au. In all 
three runs, the maximum eccentrici ty of the planetesi mals 
was set to - 0.03 , motivated by iLohne et"^ d2012h and 
ISchiippler et al. (|2014|) who showed a preference for low dy¬ 
namical excitations in debris discs around late-type stars. For the 
grain material, we chose a homogeneous mixture with volume 
filled in equal parts by astronomical silicate, amorphous carbon, 
and vacuum, denoted as Ml (Table|2]i. Other assumptions for the 


material composition, Cmax, and ArpB are probed in the Sects. l4~2l 
l431andl4Al 


4.1.1. Impact on thermal emission 


Figure |2] shows SEDs and normal optical thicknesses, t±, for 
SWO, SW50, and SW300. The vertical height of each SED was 
fitted by searching for the timestep in the ACE simulations where 
the dust mass has an appropriate value to reproduce the observed 
level of thermal emission. For the SEDs, we determined/A, 
which is the sum of the squares of the deviations of our models 
from N individual photometric points, divided by N. To this end, 
we considered the N -20 SED points for 4 > 10 pm (Table[T]i. 
We stress that the models are not fitted to the data through vari¬ 
ations of the initial disc parameters. The;^fJgp/A metric merely 
serves for a better comparison of the SED models. All models 
show good fidelity with the photometry between 160 - 850 pm 
but markedly underestimate or tend to underestimate the 1.3 mm 
data. At shorter wavelengths, SW50 provides a good match to 
the data, whereas SWO significantly underestimates the mid-IR 
points and SW300 the 70 pm point. 

In the planetesimal zone of run SWO, the amount of par¬ 
ticles continuously increases with decreasing grain size apart 
from an underabundance around s - 0.4 pm (Fig.|2}l, r — 39 au). 
Although the smallest grains in the range of a few tens of 
microns are the most frequent, they barely affect the SED as 
their thermal emission is negligible in the IR. Thus, the effec¬ 
tive minimum grain size is given by the maximum of the size 
distribution at around 10 pm. This large grain size is charac¬ 
teristic fo r discs with low dynamical excitation of the plan¬ 
etesimals (Thebault & Wul2008l : lLohne et al.l2012HKrivov et al.l 


l2013l : [^htir)pler et alJ2014 ). The size disPibutions in SW50 and 

SW300 show maxima near 0.07 pm and 0.5 pm, respectively 
(Fig. |2^,f). These values are close to the blowout sizes, Sbiow 
(Table |2]i. Most of the smaller grains move on unbound orbits 
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SWO SW50 SW300 




radial distance, r [au] 


radial distance, r [au] 


radial distance, r [au] 


Fig. 2. Disc models of narrow planetesimal belts centred at 40 au with different assumptions for the stellar wind (S W) strength (from 
left to right M* = 0,50,300 Mf). (a)-(c): SEDs. Crosses show photometric data and squares star-subtracted photometric data. The 
dashed lines depict the stellar photosphere, the dash-dotted lines the disc emission, and the solid lines the star-i-disc emission. The 
Tsed/^ were obtained using the N -20 photometric points in the plotted range, (d) - (f): normal optical thickness as a function 
of grain size at different distances from the star. The solid curves (r = 39 au) show the particle distribution within the parent belt. 
Shaded regions depict the contribution of grains on unbound orbits at r = 39 au. fg) - (i): normal optical thickness as a function of 
distance from the star. Shaded regions mark the contributions from different size ranges. 


and leave the system on short timescales. Towards larger grains, 
the size distributions in the parent belt are nearly flat until they 
become significantly steeper at sizes between 10 and 100 p.m. 
This behaviour is caused by the strong transport due to stel¬ 
lar wind drag, which favours smaller parti cles and quickly re¬ 
moves them from the planetesimal ring (e.g . .IVitense et al.l201p: 
iReidemeister et 'Zll2()l ll IWvatt et al.ir201 it iKrivov et al.ll2013h. 


As a result, the normal optical thickness is radially uniform for 
r < rpB in SW50 and SW300 (Fig.|2h,i). 


Direct stellar radiation pressure, aided by stellar wind pres¬ 
sure in SW50 and SW300, pushes submicron-sized, barely 
bound grains to eccentric orbits resulting in extended dust 
halos. However, both because of low dy namical excitation 
dThebault & Wul l2008h and fast transport (IStrubbe & Chian ^ 
l2006l) . the halos are tenuous and exhibit radial profiles with outer 
slopes that are significan tly steeper than the -1.5 predicted for 
“cano nical” debris discs (iKrivov et al.ll200^ IStrubbe & Chiangl 
I 2 OO 6 I 1 . 


For all three runs, the surface mass density of solids, S, is ris¬ 
ing with distance up to rpB and shows a bow-like shape in a log- 
log S - r diagram. Averaged over r = 1 - 30 au, we measured 


S oc for SWO, and S oc for SW50 and SW300 . The lat¬ 
ter is close to derived bv iMacGregor et ^ d2013h . The dust 
mass (objects with radii up to i = 1 r nm) is about 2 x 10~^ Mm , 
consistent with previous estimates dAugereau & Beusl l2006ll . 
The three models have a fractional luminosity of 4 - 5 x 10 
wh ich is not far fro m 6 x 10 obtained by a modified blackbody 
fit dT.iii etal.ll2004 . 

Figure [3] shows the surface brightness profile extracted from 
the ALMA 1.3 mm image. The profile increases with distance 
up to » 30 au and drops steeply beyond. Since long-wavelength 
observations reveal the spatial distribution of large particles that 
have similar orbital elements as the planetesimals, this break in¬ 
dicates the outer edge of a dust-producing planetesimal zone. 
For modelling of thi s profile, we assu med a disc inclination of 
89.5° from face-on dKrist et al.l l2005h and convolved the syn¬ 
thetic 1.3 mm images for SWO, SW50, and SW300 with the 
ALMA reconstructed dirty beam as produced by CASA after 
proper rotation. From these images, we extracted the radial pro¬ 
files by integrating over a 0.1" central strip along the disc ma¬ 
jor axis. We allow for vertical scaling of the gathered synthetic 
profiles. The scaling factors account for an imperfect reproduc- 
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Fig. 3. ALMA radial surface brightness profile at 1.3 mm (dots) 
vs. models (lines). Synthetic profiles are multiplied by the fac¬ 
tors given in the legend to compensate imperfections of the mod¬ 
elling procedure (explained in the text). 


tion of the total flux density at 1.3 mm, visible by comparing the 
models with the observed SED. Such imperfections stem, e.g., 
from uncertainties in the absolute flux calibration or in the op¬ 
tical properties of the dust grains (iLohne et alJl2012h . Further, 
they derive from a real inability of our models to reproduce the 
ALMA flux density along with the other photometric data. We 
surmise that the latter reason mainly causes the models to dis¬ 
tinctly underestimate the ALMA flux density and necessitates 
high scaling factors (Fig. |2]l. In addition, the dust model used 
possibly underestimates the dust emissivity at mm wavelengths. 
Figure [3] clearly shows that the surface brightness peak moves to 
smaller projected distances with increasing stellar wind strength, 
reflecting the enhanced inward transport of particles. 

4.1.2. Impact on scattered light 

We calculated the surface brightness of scattered light for our 
models along the projected distance, b, in the disc midplane by 

cn^ T r r ^■5^2sca5ll(0)n(r, S) 

S(b)^L, J J -- dsdl, (2) 

where L* is the luminosity of the star, n(r, s)ds the dust num¬ 
ber density of grains with radii in [s, s + ds] at a distance r, and 
I = +(r^ - b^y^^ the line of sight. Qsca and S ii{(p) are the scatter¬ 
ing efficiency and the scattering phase function of the grains, re¬ 
spectively. The latter can be retrieved from the 4 x4 Mueller ma¬ 
trix S ij and is a function of the scattering angle (ft - arcsin(Z7/r). 
We used Mie theory to determine Qsca and S ii(0). The ACE sim¬ 
ulations provided the dust number density. 

To compare with the o bserved V- and //-band surface bright¬ 
ness profiles presented in iFitzgerald et al.l (l2007h . we integrated 
S (b) over the filter transmission curves of the HST Advanced 
Camera for Surveys F606W (central wavelength Tc = 0.59 |xm, 
width A/I = 0.23 pm) and the Keck NIRC2 camera (Ac - 
1.63 pm, A/1 = 0.30 pm). The disc surface brightnesses were 
expressed relative to the star for which we found 8.55 mag in 
V-band and 4.84 mag in //-band with our photosphere model. 
To account for the extended halo of small particles, the flux was 
integrated up to stellocentric distances of r = 150 au for each 
line of sight. We also compared with the radial p rofiles extracted 
from the high-resolution HST/STIS image by ISchneider et alJ 


(l2014l) . In their Fig. 39, the surface brightness is given for both 
sides of the disc along with their exponentional fits. We averaged 
both fits and converted the flux from mJy/arcsec^ to mag/arcsec^ 
using the zero-point magnitude flux of 3671 Jy and then sub- 
ttacted a s tellar V-band br i ghtnes s of 8.837 mag according to 
Table 1 in ISchneider et al.l (12014h . The resulti ng V-band pro¬ 
file wa s found to be in good agreement with the lFitzgerald et al.l 
(l2007l) data. 

Figure |4] depicts the modelled and observed V- and //-band 
profiles as well as the disc colour, V - H. All models show 
significant deviations from the surface brightness data in both 
bands. Although roughly reproducing the observed slopes for 
b <40 au, the synthetic profiles are steeper beyond, which is 
the halo zone, but get shallower with increasing stellar wind 
strength there. In terms ofV-H, the SWO and SW50 discs ap¬ 
pear slightly blueish (V - H < 0), whereas the SW300 disc is red 
(V - H > 0), which completely disagrees with the observations. 
These effects are mainly due to the distribution and the amount 
of small grains in the disc. Increasing the stellar wind strength 
pushes more small particles in the halo zone (Fig. |2]i, and helps 
to enhance the scattered light flux beyond 40 au to approach the 
observed level. However, also the blowout size increases and the 
whole disc becomes more and more populated by larger parti¬ 
cles that are stronger scatterers at longer wavelengths. Hence, 
the disc colour V-H switches to red for too strong winds, clearly 
conflicting with the observational data. 

Next, we considered the degree of linear polarisation, given 
by 


P(b)^ 


f f p((f>) 5ii(<^) TTi^gsca n(r, i) di dl 
f f Sn(<P) 2sca n(r, s) ds dl 


(3) 


where p((p) is the linear polarisation of a single particle. If P(b) is 
positive (negative), the scattered light is partially polarised per¬ 
pendicularly (parallel) to the scattering plane, i.e. the plane con¬ 
taining the star, the dust grains, and the observer. 

As shown in Fig. |5] the observed monotonous increase of 
the degree of polarisation up to 80 au is roughly reproduced by 
SW50 only. In run SW300, the halo region is dominated by par¬ 
ticles of 0.5 pm in size which are about one order of magni¬ 
tude larger than in SW50. Their polarisation p(^) strongly oscil¬ 
lates with the scattering angle. Accordingly, the polarisation in¬ 
tegrated over a range of scattering angles along the line of sights 
tends to average to small values and P{b) does not rise with dis¬ 
tance. For SWO, P(b) increases up to the planetesimal belt as 
observed, but then decreases beyond. 

Although SW50 fits the measured degree of polarisation 
best, the data are markedly overestimated by this model for 
b <40 au. This is ca used by particl e s drag ged from the parent 
belt to the inner zone. iGraham et akl (12007 ) and iFitzgerald et akl 
(l2007h concluded that the inner disc region has to be relatively 
free of scattering grains, resulting in a low normal optical depth 
to scattering, - Jns^ Qsc3.n(r, s)ds. In V-band, SW50 

gives ry^(r < rpB)/T^'‘(r = rpe) > 0.1, which is a t least two or¬ 
ders o f mag nitude greater than the limit found by iGraham et al.l 
(l2007h and iFitzgerald et alJ (l2007l) . However, owing to strong 
stellar winds, the inward transport of a significant amount of 
dust is a natural outcome in the dynamical evolution of the outer 
planetesimal belt. Even if planets exist at r < rpp, small particles 
come through, as only slower, bigger particles can be efficiently 
scatt ered by planets or get trapped in mean-motion resonances 
(e.g.. lReidemeister et al.l20rillshannon et al.ll2015l) . It is impor¬ 
tant to note that the deviations of the model from the data may 
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Fig. 4. Scattered light profiles of the AU Mic disc, (a) - (c): blue and red lines show synthetic midplane surface brightness prohles 
in V- and //-band relative to the stellar brightness for the disc models SWO, SW50, and SW300 (Sect. 14.11) . Curves labelled with 
Mie in the legend were calculated with scattering phase functions obtained by means of Mie theory. The blue dou ble-dotted line 
in pan el (b) depicts the V-band prohle generated with the best-ht Henyey-Greenst ein (HG) phase function model of iGraham et al.l 
(l2007h . Dots represent measurements with Icr errors from iFitzgerald et^ (l2007h . abbreviated as F07, averaged over the SE and 
N W wings (circles for V- band, squares for //-band). The black dashed line shows the averaged V-band ht to both wings obtained 
bv ISchneider et al.l (l2014l) . abbreviated as S14. (d) - (f): observed and modelled colour prohles, V — H, obtained with the F07 data 
and the Mie prohles shown in panels (a) -(c). 



Fig. 5. Degree of polarisation as a function of projected dis¬ 
tance. Squares and tria ngles are measurem ents for the SE and 
NW ansae, respectively (iGraham et al.l200^ . The double-dotted 
green line wa s obtai ned with the Henyey-Greenstein model of 
IGraham et al.l (l2007h . whereas the other lines are results from 
Mie theory. 


be caused by the assumption of spherical dust grains. Real parti¬ 
cle shapes signihcantly deviate from spheres because the bodies 
underwent a rich accretional and collisional history. While effec¬ 
tive medium Mie spheres provide a good approximation to the 
total scattering cross sections of irregular particles, they cannot 
reprod uce the scattering phase function and polarisation of such 
grains (IShen et al.l2008[[20091) . Hence, Mie theory fails in repro¬ 
ducing the scatt ered light observations as it was rece ntly shown 
for HR 4796 A (iMilli et ahllMl iPerrin et alJ l2015ll. In a study 
of the HD 181327 disc. IstSc et al. ( 201 dll also found an em¬ 
pirical scattering phase function, which is not fully reproducible 
with Mie theory, indicating strongly forward scattering grains. 
Thus, if the AU Mic disc is populated by similar particles, our 
models would poorly ht the scattered light measurements by us¬ 
ing Mie theory. 


In an additional test, we aim at illustrating the effect when 
deviating fr om the Mie sphere d ust model. To this end, we 
adapted the IGraham et alJ (l2007l) best-ht Henyey-Greenstein 
(HG) model for 5ii(0) and p(0), and re-computed the V-band 
surface brightness prohle and the degree of polarisation for 
SW50 (double-dotted curves in Eig. |4jr and Eig. |5]). Note that 
5 11 (0) and piip) of the HG model are independent from grain 
size. Thus, all disc objects have the same scattering properties. 
The new synthetic U-band prohle hts better the absolute val¬ 
ues of the surface brightness Hux for h < 40 au, whereas the de- 
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gree of polarisation rises too steeply in that range. Reasons for 
this behaviour can be identihed by comparing 5ii(0) and p{(p) 
predicted by Mie theory and the HG model (Fig. 1^. The HG 
phase function is distinctly broader than for Mie spheres with 
s > 0.5 pm. That is why bigger particles can scatter more light 
towards the observer that leads to the V-band flux enhancement 
for h < 40 au as shown in Fig. |4]r. Beyond 40 au, only small par¬ 
ticles are present and the HG profile resembles the one of the 
Mie spheres. 

However, the HG model shows strong deviations from the 
observed degree of polarisation, attributable to the form of the 
polarisation curve as a function of scattering angle (Fig. ig. The 
maximum of the polarisation occurs at (^max = 90°, so that the 
main contribution to the polarised light comes from particles 
with distances from the star nearly equal to their projected dis¬ 
tances (r ss b). The of the Mie sphere polarisation func¬ 
tions strongly vary with s. T his agrees with th e trends found for 
irregularly-shaped particles (IShen et al.ll2009h . With such shifts 
of ipmax the polarised light at a certain b mainly originates from 
particles at r > b. Furthermore, irregular grains smaller than 
0.1 pm tend to be weaker polarisers than Mie spheres in the same 
size range. Both may help to generate a shallower increase of the 
degree of polarisation with projected distance. Thus, the consid¬ 
eration of irregular particles in the scattered light analysis seems 
to be an attractive possibility to mitigate the deviations between 
the observed and modelled degree of polarisation. 

We summarise that all our models have a large amount of 
small grains within the parent belt in contrast to what was found 
in previous studies and none of our models provides good fits 
to the scattered light data. There are two ways to interpret this 
result: (1) Our models may predict the actual dust distribution 
in the disc well but Mie theory is not valid to simulate the scat¬ 
tering properties of the grains. Then, our analysis show a prefer¬ 
ence towards a scattered light model where the grains are weaker 
polarisers than Mie spheres with polarisation maxima at scatter- 



0 30 60 90 120 150 180 


scattering angle, ^ [deg] 


Fig. 6. Normalised scattering phase function (top) and polar¬ 
isation (bottom) vs. scatte ring angle in V-ban d. Solid black 
lines have been inferred bv iGraham et al.l (l2007h via a Henyey- 
Greenstein model. The coloured lines are for Mie spheres with 
radii s, made of material Ml. 


ing angles different from 90°. (2) Mie theory is an appropriate 
technique for the scattered light analysis but there are shortcom¬ 
ings with respect to the dust density derived in the collisional 
modelling. The truth may also be a combination of (1) and (2). 
A consideration of point (1) goes along with the simulation of 
the emission from irregular particles which is beyond the scope 
of this work. Thus in the following, we consider only point (2) 
in more detail. Since the scenario SW50 is found to be in best 
agreement with the data amongst the three stellar wind strengths 
tested, we consider SW50 a reference model, which is the start¬ 
ing point for further individual parameter variations influencing 
the dust distribution. 


4.2. Material composition 


The chemical composition of the AU Mic disc might be simi¬ 
lar to discs around other BPMG stars because they have formed 
in the same molecular cloud from the same material and are 
coeval. However, the disc objects are heavily processed during 
the debris disc and preceding protoplanetary disc phases. Thus, 
their composition may depend on the mass and luminosity of 
the host star. Circumstellar dust compositions for s ome BMPG 
memb ers h ave already be e n ana lysed previously. ISmith et^ 
(l2009h and IChurcher et alJ (1201 Ih considered grains with a sil¬ 
icate core and a mantle of organic refractories for SED mod¬ 
elling in the rj Tel and HD 191089 systems. They pointed 
out that particles with small silicate cores and a porosity of 
20% (rj Tel) and 60% (HD 191089) provide good fits to the 
observed excesses. A similar model that incorporates porous 
core-mantle silicate grains, additionally covered by ice in the 
outer region of the disc, also reproduces the B Pic observa¬ 
tions dPa ntin et al. 1997t ETife Greenberdl 19981 : lAugereau et al.l 
l2001h . Lebreton et al.l ( ^12ll tested various chemical composi¬ 
tions for the debris around HD 181327. They found the SED to 
be mostly consistent with grains consisting of amorphous sili¬ 
cate and carbonaceous material with a dominant fraction of ice. 
In addition, the particles have to have a porosity of 65%. 

The material composition affects the dust distribution and 
the SED by changing /? an d the dust temperature, (e.g., 
iKrrchschlager & Wolf 120131 for the influence of porosity on p 
and Td). We now depart from the mixture Ml, used in our ref¬ 
erence run SW50, to assess whether other materials can be vi¬ 
able for the AU Mic system. We repeated the simulation with 
the compositions M2, M3, and M4 explained in Table |2] and 
investigated modifications to observables (Eig.|7|i- 

The material M2 represents compact particles, composed of 
silicate and carbon. Eor these, we found a shift of the SED to¬ 
wards shorter wavelengths that deteriorates the fit in the submil¬ 
limetre to radio range. In M3 and M4, we included water ice 
additionally. All constituents are present in equal volume frac¬ 
tions in M3, whereas the porosity and ice fraction are increased 
in M4. M3 is in fairly good agreement with the observed pho¬ 
tometry, only tending to underestimate the near-IR excess. The 
icier and more porous M4 produces an SED whose peak posi¬ 
tion remains rather unchanged but increases and decreases more 
steeply. This narrower SED distinctly disagrees with the near-IR 
and submillimetre measurements. 

No large differences in terms of the degree of polarisation are 
noticeable for Ml, M2, and M3, whereas M4 markedly underes¬ 
timates the data. The latter gives evidence that the material must 
not be too icy and porous. Note that modelling of the degree of 
polarisation is possibly affected by limitations of the Mie theory 
(see Sect. 14.1.21) that does not allow us to find strong constraints 
for the upper limit of the ice and vacuum content. 
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wavelength [pm] 



Fig. 7. Influence of the dust composition on the SED (left) and the degree of polarisation (right). The material names Ml - M4 are 
explained in Table For better comparison of the SED models, /A values are given. 


We noticed little changes in the radial profile at 1.3 mm for 
all materials considered. This is caused by two characteristics 
of large grains, dominating the long-wavelength emission. First, 
large grains are weakly affected by stellar radiation and stellar 
winds and their spatial distribution is equal for different materi¬ 
als. Second, the dust temperature approaches the blackbody tem¬ 
perature with increasing parti cle size at a given distance (e.g.. 
Fig. 5 of iPawellek et al.ll2014l) . Therefore, Td is nearly indepen¬ 
dent of the chemical composition for large grain sizes. 

Altogether, the SED and the degree of polarisation seem to 
be consistent with materials having a small fraction of ice and a 
moderate porosity. 

4.3. Dynamical planetesimal excitation 

The maximum eccentricity of the planetesimals, Cmax, represents 
the disc’s dynamical excitation since this parameter affects the 
impact velocities of the objects, and therefore, the amount of 
small grains produced in collisions. In addition to the refer¬ 
ence model with e^nx - 0.03, we consider e^nx - 0.01, a low 
value corres ponding to the expec ted maximum level of pre¬ 
stirred discs (iMatthews et al.ll20l4 and references therein), and 
Smax = 0.1, a high value representing the dynamical excitation in 



Fig. 8 . SED models for different values of the maximum eccen¬ 
tricity, Cinax, of the planetesimals orbits. 


the Edgeworth-Kuiperbelt (e.g.. lVitense et alJ2010h . The result¬ 
ing SEDs are shown in Fig. [8] 

Increasing e^ax favours the production of small particles 
and shifts the SED maximum to shorter wavelengths. Thus, 
the model with e^ax = 0.1 well agrees with the excess data at 
A <70 pm, but completely fails to reproduce the ones at longer 
wavelengths. In contrast, e^ax - 0.01 provides a good match 
to the photometry beyond 100 pm, but strongly deviates from 
the measurements at d < 70 pm. Improvements may be possi¬ 
ble by additional variations of other parameters. For instance, 
shifting the planetesimal belt inwards (outwards) increases (de¬ 
creases) the dust temperature and helps to And a better SED 
fit for Cmax = 0.01 (Cmax = 0.1). Howevef, any significant rpB- 
shift would result in stronger deviations from the ALMA 1.3 mm 
profile, which constrains the location of a narrow planetesimal 
belt at around 40 au. Another possibility would be broadening 
the SED towards shorter wavelengths for e^ax = 0.01. This can 
be achieved by stronger stellar winds as depicted in Fig. |2] 
However, sponger winds also lead to an increase of the blowout 
size, and therefore, it is more difficult to reproduce the observed 
blue disc colour (Sect. l4~TT2l i. Furthermore, other dust materials 
may have a positive effect on the SEDs. None of the materials 
investigated in Sect. I4.2l spengthens the submillimePe flux den¬ 
sity, as required to improve the SED for e^ax = 0.1. Only one 
composition - the non-porous M2 material - predicts a stronger 
mid-IR emission from which the low-Cmax model would benefit, 
but this material agrees less with the scattered light measure¬ 
ments. In summary, the reference run with Cmax = 0.03 matches 
the entire SED the best. Although values higher or lower than 
0.03 seem to be unlikely, they cannot be ruled out completely if 
combined with other parameters, which was not explored here. 

An observational constraint on the dynamical excitation can 
be put by the vertical thickness of the disc through the equi¬ 
librium condition i - e/2, where i and e are the orbital in¬ 
clinations and eccentricities of the disc particles. Thus, mea¬ 
suring the disc’s opening angle at long wavelengths hints at 
the maximum eccentricity of the planetesimals. However, the 
AU Mi c disc is vertically unresolved in the SMA and ALMA 
images (IWilner et al.l2013;lMacGregor et al.ll2013h . which gives 
only upper limits of the disc extent in this direction. We 
compared synthetic vertical profiles of our e^ax = 0.1 model 
with the ones exPacted from the ALMA image. The modelled 
profiles are marginally consistent with the ALMA data, con- 
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straining the disc’s semi-openin g ang k to /max ^ gmp/2 = 0 .05. 
Interestingly, iKrist et aTI (l2005l) and iMetchev et alJ (l2005h de¬ 
rived imax = 0.04 - 0.07 from visible and n ear-IR images. 
Readin g off the vertical scale height from the ISchneider et al.l 
(|2014|) high resolution STIS image gives /max = 0.03. Thus, the 
observations yield a maximum eccentricity Cmax = 0.06 - 0.14 
that is excluded by our modelling. Resolving this contradic¬ 
tion may be possible if the vertical disc thickness seen in short- 
wavelength images does not neces sarily points tow ards the dy¬ 
namical state of the nlanetesimals. iThebauld (l2009h proposed a 
possibility for the AU Mic disc to develop a vertical thickness 
close to the observed value in the visible and near-IR although 
the disc is dynamically cold. In that model, the smallest grains 
have high in-plane velocities due to radiation and stellar wind 
pressure, and these velocities are partially converted into ver¬ 
tical ones in collisions with other particles. As a consequence, 
the disc becomes naturally thicker even without the gravitational 
perturbation of large, massive embedded bodies. In equilibrium, 
the disc has a large vertical dispersion of the smallest grains, 
whereas large particles remain close to the midplane. Thus, a low 
^max, as preferred in our collisional modelling, can be consistent 
with the observed vertical disc thicknesses at short wavelengths. 




4.4. Planetesimal belt width 
4.4.1. 17au-wide belt 

As noted in Sect. 14.1.11 the outer edge of the planetesimal belt 
is well constrained by the ALMA profile. However, the edge-on 
orientation of the disc and the limited information on the inner¬ 
most regions affected by the bright central source of unknown 
nature make the position of the inner edge of the disc, and so 
the disc extent, quite uncertain. In another test, we increased the 
planetesimal belt width from ArpB = 5 au (reference model) to 
ArpB = 17 au by shifting the inner edge to 25.5 au while keep¬ 
ing the outer edge at 42.5 au. The initial surface density was as¬ 
sumed to be constant over the belt width. Figure |9] compares the 
17 au-wide belt scenario with the reference model. 

The broader planetesimal belt has nearly the same SED as 
in the reference run and does not increase the warm emission 
although the dust production zone is more extended towards the 
star. The 1.3 mm profile of the 17 au-wide belt peaks closer to 
the star than for the reference model. Thus, it tends to under¬ 
estimate the observed surface brightness farther out. This short¬ 
coming does not hint at an implausibility of a broader planetes¬ 
imal belt model since it can easily be compensated by shifting 
the outer edge from 42.5 au outwards. Although the broader belt 
scenario matches better the V-band profile than the reference run 



Fig. 9. Influence of the planetesimal belt width, Arpp, on (a) the SED, (b) the 1.3 mm radial profile, (c) the V- and //-band profiles, 
and (d) the degree of polarisation. The filled patterns in panel (b) depict the location and the radial width of the planetesimal 
belts (black-hatched 5 au, red-shaded 17 au). The profiles were vertically scaled by 1.7 (ArpB = 5 au) and 1.5 (ArpB = 17 au) to 
compensate the modelling imperfections mentioned in Sect. 14.1.11 


13 






























































Ch. Schiippler et al.: Collisional modelling of the AU Microscopii debris disc 




Fig. 10. Left, time evolution of the surface density in an extended planetesimal disc (ArpB ~ 45 au). The wavy structures of the 
curves are numerical artefacts. For the sake of comparison, the straight line shows the Z-slope derived bv iMacGregor et ^ (l2013h . 
Right: time evolution of the 1.3 mm profile along with the ALMA data. Note that the 0.01 Myr surface brightness profile was scaled 
by a factor of 0.05 for better visibility. 


does, the degree of polarisation shows stronger deviations which 
may be attributed to the light scattering model (Sect. l4~L2l i. 

4.4.2. 44au-wide belt 

We now address the possibility that the AU Mic disc was pop¬ 
ulated by planetesimals over a wide range of distances. In such 
a disc, the innermost planetesimals have the shortest collisional 
lifetimes and are depleted within a few Myr, resulting in a sur¬ 
face density increasing outwards (iKennedv & Wvat3l2010h . Due 
to the inside-out erosion of objects, an extended planetesimal 
disc may look like a narrow belt in long-wavelength observa¬ 
tions, even if there are no planets within the disc that have swept 
up the inner region. In their study of the collisional evolut ion of 
a fiducial, lOOau-wide disc. iThebault & Augereaul (l2007h con¬ 
firmed that the surface brightness peaks at (sub)millimetre wave¬ 
lengths close to the outer edge of the planetesimal distribution 
after several Myr. Thus, the extended-disc scenario is an auspi¬ 
cious alternative for explaining the rising ALMA profile up to 
projected distances of 30 au in the AU Mic system. The signif¬ 
icant fall beyond 30 au may mark the point beyond which plan¬ 
etesimals have not formed at all or where the disc was truncated. 

To assess the feasibility of this scenario, we started an ACE 
simulation with a broad radial distribution of planetesimals from 
1 to 45 au. Here we already accounted for a slight shift of the 
outer edge from the reference value 42.5 au to 45 au, as sug¬ 
gested in Sect. 14.4.11 The run was initialised with a solid sur¬ 
face density given by the Minimum Mass Solar Nebula (MMSN) 
model, We considered bodies up to 100m in ra¬ 

dius and assumed pre-stirred planetesimals with Cmax = 0.03 
(Sect.l43Ti. 

Figure [10] depicts the evolution of the surface density and 
the 1.3 mm profile. The surface density quickly dropped close to 
the star due to particle erosion and the initially negative slope 
(-1.5) switched to a positive slope in a few ten thousand years. 
Note that the particles’ collisional lifetimes were shorter than 
their orbital periods at the beginning of the simulation due to 
the large particle density at that phase. Thus, two consecutive 
collisions of one object are correlated and the orbit averaging 
method, implemented in ACE, is a crude approximation to deter¬ 
mine the phase space density. Therefore, the results for the first 


one hundred thousand years of simulation time should be treated 
with caution. 

Owing to the progressive inside-out depletion of the plan¬ 
etesimals, the maximum of the 1.3 mm profile moved towards 
larger distances as time elapses. After 3 Myr, we noticed that 
the profile peak reached a maximum distance of about 30 au, 
which is in good agreement with the ALMA data. The absolute 
flux level of the observed dust emission was reproduced after 
29 Myr. The final run fits the SED, the scattered light profiles, 
and the degree of polarisation with nearly the same fidelity as 
the 17 au-wide belt model (Sect.|4j4Tj. 

Note that the collisional evolution of the disc was artificially 
shortened as the largest bodies in the disc had been 100 m in size. 
If the disc is contains larger bodies, the inside-out planetesimal 
depletion is slower and the disc needs m ore time to reach the 
collisional equilibrium (iLohne et al.ll2008l) . As a result, it takes 
longer than 3 Myr to reach the observed shape of the 1.3 mm pro¬ 
file. This way, scaling the model to the proper age of the system 
is possible. However, we refrained from further simulations, be¬ 
cause of the numerical complexity of the extended planetesimal 
disc simulatiorOl The run presented does illustrate the feasibility 
of an extended planetesimal disc scenario sufficiently well. 

5. Unresolved central emission 

In Sect. |4| we did not care about the nature of the inner unre¬ 
solved component and assumed its contribution to the flux den¬ 
sities at all wavelengths to be negligible. However, this may not 
be true if the unresolved central emission originates from an in¬ 
ner dust ring. The dust there would be much warmer than in 
the outer disc and would emit significantly at short wavelengths, 
while adding little flux in the mm range. We now consider the 
possibility of an additional inner planetesimal ring in the AU Mic 
system. If such a ring could exist, we should also investigate its 
consequences for the modelling of the resolved outer disc, dis¬ 
cussed in Sect. |4| 

A putative inner dust ring has to be 6 times brighter than 
the stellar photosphere at 1.3 mm (iMacGregor et al.ll2013L and 
Sect. 12.21 of this paper). Simultaneously, the ring has to emit 

' It took about 150 CPU days to evolve the disc over 29 Myr. 
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weakly in the mid-IR in order not to contradict the measured 
excess there. Since an ALMA resolution of about 0.6” («; 6 au) 
was achieved, the ring should be located at distances < 3 au to be 
compatible with remaining unresolved. We placed an inner plan- 
etesimal ring at rpB,in = 3 au, having a full width ArpB,in =0.4 au, 
and assumed Cmax = 0.03 and 50 times the solar wind strength. 
We made sure that the ring emission did not exceed the excess 
shortward of 24 pm and then investigated its emission at 1.3 mm. 
Figure [TT] shows the results. 



Fig. 11. SEDs of the inner planetesimal rings centred at 
^PB.in = 3au (width ArpB.in = 0.4au) for the compositions Ml, 
M3, and M5 (Table |2ll. Squares illustrate the photosphere sub¬ 
tracted photometry. The black cross in the right lower corner 
shows the emission of the unresolved star subtracted central 
component at 1.3 mm with Icr uncertainty. The dashed line rep¬ 
resents the M3 outer disc in Fig.|7]in combination with the inner 
M5 disc. 

First, we assumed two inner rings with Ml and M3 - ma¬ 
terials we found suitable for the outer component (Sect. lO) . 
Consistency with the short-wavelength excess was achieved with 
a maximum dust mass of 5 x 10“^ Me. The rings are not bright 
enough to completely reproduce the central emission. Although 
M3 enhances the SED maximum, the long-wavelength slope is 
steeper and levels to nearly the same flux density at 1.3 mm as 
for Ml. In another attempt, we considered highly porous grains 
(M5). This is motivated by the possibility that icy grains are 
produced beyond the ice line and then drift towards the star by 
drag forces. Accordingly, the particles get heated and their sub¬ 
limated ice fractions are replaced by vacuum inclusions, which 
increase the porosity. Assuming an ice subl imation tempera¬ 
ture Tsubi = lOOK (e.g.. lKobavashi et al.ll2008l) yields an iceline 
at Hce = (Tsubi/277 K)“^ au 2 au, which is close to 

^PB.in- The M5 material emits slightly more than Ml and M3 in 
the mm regime, but reproduces only about 60% of the observed 
central flux density. This constrains the maximum fraction of 
the central emission that can be explained by pure dust emis¬ 
sion. Thus, the central emission must originate only or in parts 
from another source. Note that the star itself can contribute to 
the central emission, for instance due to an active stellar corona 
jCranmer et al.l 12013 . Furthermore, stars may emit more flux 
at long wavelengths than predicted by the purely photospheric 
models for their spectral class. ALMA observations of the G- 
and K-type stars a Cen A and B show a flux enhancement by 
20-30% in band 7 and almost doubling of the flux in band 3 


(iLiseau et al.ll20L5h . Provided such a trend continues to an M 
dwarf like AU Mic, we do expect a significant stellar contribu¬ 
tion to the observed central emission. 

To get the total SED, the inner ring should be combined with 
an appropriate model for the outer disc. In order to avoid con¬ 
tradictions with the data, the outer disc has to emit less flux in 
the mid-IR where the SED of the inner ring peaks. This is well 
achieved by the M3 outer disc (Sect. lO l. and indeed, in com¬ 
bination with the M5 inner ring the overall SED is well repro¬ 
duced (Fig.[TTl dashed line). Furthermore, a combination of the 
icy material in the outer ring and iceless, highly porous material 
in the inner ring is reasonable, as explained above. We adopted 
the outer disc SED from the model described in Sect. 14.21 with- 
out any adjustments (e.g., vertical scaling) that one possibly ex¬ 
pects due to the presence of the inner ring. Thus, even if the in¬ 
ner emission stems from dust, it seems to barely affect the outer 
disc modelling. The^gg^/A increases from 2.7 for the M3 outer 
disc to 3.3 for the combined M5 inner and M3 outer disc. The 
larger /A is mostly due to the deviation of the inner-i-outer 
disc model from the Spitzer/MlPS 24 pm point. However, the 
overestimated flux density at 24 pm may be compensated by a 
modification of the dust composition. As seen in Fig. [T] adding 
ice reduces the flux in the mid-IR, while leaving it nearly un¬ 
changed at longer wavelengths. This promises improvements if 
a material icier than M3 for the outer disc is used. We highlight 
the small shift of the blowout size when going from iceless to 
icy grains (Tab. |2]i. Accordingly, the abundance of small grains 
will not change significantly, still reproducing the scattered light 
data (Sect. l4~T!2l l. 

We stress that our two-component model does not account 
for the collisional interaction between inner and outer compo¬ 
nent which is mainly caused by two effects: (i) collisions of in¬ 
ner disc particles with inwards dragged grains from the outer 
disc and (ii) collisions between unbound particles produced in 
the inner ring with outer disc particles. By comparing the nor¬ 
mal optical thicknesses of inner and outer component, including 
bound and unbound grains, we found Tj^_out/T±,in ~ 0.1 at 3 au 
and Tj^,in/T±,out ~ 0.001 at 40 au. Thus, collisional interaction be¬ 
tween both is weak, justifying our modelling strategy. 

6. Comparison with parametric modeiiing 

Collisional modelling allows us to explore only a small number 
of parameter combinations. Although the models found repro¬ 
duce the data satisfactorily, they are not unique, and better so¬ 
lutions may exist. In order to And an independent confirmation 
of the collisional modelling results, we additionally analysed the 
data with a multidimensional fitting approach. In that approach, 
the dust number density was assumed to be a combination of two 
independent power laws, representing the size and radial distri¬ 
bution of the grains. This model does not consider the physical 
mechanisms of the dust production and evolution. However, it is 
well suited for an exploration of a wide parameter space to find 
the most likely solutions that we compared with the collisional 
modelling results. 

In the power-law model, the surface density follows r“ be¬ 
tween the inner and outer cut-off radii, R^n and Rout, while the 
grain size distribution follows between the lower and upper 
cut-off grain sizes, Smin and Smax- The influence of Smax is negli¬ 
gible for reasonably steep grain size distributions with y < -3. 
We fixed i^ax to 2 mm. The dust mass Md is an additional free 
parameter. The disc inclination was fixed to 0 = 90°. For the dust 
composition, we assumed co mpact, spheric al grains composed 
of pure astronomical silicate (lDrainell2003h . This makes a total 
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Table 4. Power-law fitting results. 


Parameter 

Range explored 

# values 

Distr. 

Best-fit [3cr] 

+ inner 

Rm [au] 

2-60 

187 

temp 

3.1* [2.0*-24.8] 

14.5 

Rom [au] 

3-100 

155 

temp 

41.9 [37.3-45.5] 

41.9 

a 

-5-5 

101 

lin 

2.8 [1.8-4.8] 

3.5 

■Smin [pm] 

0.2 - 20 

494 

log 

0.2* [0.2* - 0.5] 

1.8 

iiuM [pm] 

2000 

1 

fixed 



7 

-5--3 

21 

lin 

-3.3 [-3.4--3.3] 

-3.3 

6 

90° 

1 

fixed 



Mi [Me] 



cont 

4.7 X 10^^ 

5.7 X 10 -^ 


Notes. The distributions of the values considered in the parameter space are: temp - by equal steps in temperature of the grain size with the steepest 
radial temperature gradient, lin - linear, log - linear in the logarithm of the parameter, - fixed value (no distribution at all), cont - continuous 
(scaling of the disc mass to minimise the;y^ for given values of all other parameters). Values marked with * are at (or very close to) the boundaries 
of the parameter space explored and cannot be considered as reliable. The real values may lie outside the parameter space explored. The column 
“+ inner” shows best-fit values of the outer disc after subtracting a model of the inner, unresolved component, assuming it to be a debris disc. 
These values illustrate the uncertainties of the model parameters due to the unknown nature and uncertain parameters of this component. The mass 
values given are for grains up to i- = 1 mm. 


of eight free parameters, R\^, /?out, oi, Smin, 7, the dust com¬ 
position (sampled by only one value, namely 100% astronomi¬ 
cal silicate), and the vertical scaling of the ALMA profile. For 
the fit ting of the SEP and the ALMA prof ile, we used the SAnD 
code (lErtel et al.l2012aHL6hne et al.ll2012h which was originally 
developed for fitting simultaneously SED data and spatially re¬ 
solved Herschel data as part of the //ersc/teZ/DUNES modelling 
tool box. Note that we did not include the scattered light data for 
this power-law fitting. 

In a first approach, we fitted the model to the data ignoring 
the contribution of the central component. The best-fit model 
(TablelH reveals features we already found, at least qualitatively, 
by collisional modelling: an outward-increasing surface density 
with an index close to 2.8, a well constrained outer disc radius 
at about 40 au, a small lower grain size < 1 pm consistent with a 
moderate level of stellar winds, and a relatively shallow slope of 
the grain size distribution, y - -3.3, indicating that small parti¬ 
cles are more affected by transport than collisions. Eurthermore, 
the results are widely consistent with earlier stud ies of the disc 
dAugereau & Beu^l2006tlMacGregor et alJl2013h . 

In a second approach, we assumed the inner component to 
be an additional ring of dust. The ring was placed at 2 au from 
the star with an extent of 0.2 au and a constant surface density, 
which is just compatible with being mostly unresolved by the 
ALMA observations. The lower grain size, Smin.in, and the expo¬ 
nent, 7 in, of the size distribution for this inner component were 
explored by hand. The model was scaled to the flux density of 
the central component seen in the ALMA data. We found a very 
large lower grain size, Smin.in = 5 pm, as well as a very flat grain 
size disPibution, yin = -3, to be necessary in order to keep the 
dust cold enough (the grains large enough) not to produce too 
much excess at mid- to far-IR wavelengths. It is important to 
note that this model serves to estimate the potential contribution 
of the unknown inner component, which is probably not a debris 
disc, to the SED. However, the fitting result, even the large lower 
grain size, might be consistent with a normal debris disc. A large 
■Smin.in IS naturally expected if the grains do not experience a radi¬ 
ation pressure blowout limit as is the case for AU Mic. Then, the 
dominant or critical grain size, Sc, is the one where collisional 
and transport lifetimes are equal, which can be reached at much 
lar ger sizes. Including the effect of stellar wind drag in Eq. (6) 
of Kuchner & StarkI (120101) . and assuming an optical depth of 


10“^ at 3 au in accordance with our M5 inner ring model, yields 
a critical size Sc ~ 3 pm, which is not far from Sniin.in = 5 pm. 

In a final step, we subPacted the contribution of the inner 
ring from the SED data and re-fitted the remaining flux by the 
two-power-law model used for our first approach (Table H] last 
column). Comparing with the previous results illustrates the un¬ 
certainties of our first approach. The values are almost all within 
the 3cr-confidence interval of the first fit. This demonsPates that 
the inner component marginally influences the modelling of the 
resolved outer disc, as already found by collisional modelling. 

7. Conclusions 

We have performed in-depth collisional modelling of the 
AU Mic debris disc for the first time. A wealth of observational 
data have been considered, including the densely sampled spec- 
Pal energy disPibution, the ALMA 1.3 mm thermal emission 
profile, the scattered light profiles in V- and //-band, and the 
degree of scattered light polarisation as a function of projected 
distance in the disc midplane. Our study presents the first at¬ 
tempt to reproduce scattered light observations of a debris disc 
with collisional modelling. 

Our models provide generic radial and size distributions of 
the particles in the whole disc. However, the disc also possesses 
several subsPuctures and asymmePies like bumps and warps 
(iLiu et alJ l2004 ISchneider et al.l 1201 4l) that may be caused by 
the perturbing effects of so far undetected planets. As our mod¬ 
elling technique is limited to axisymmePic discs, the formation 
of such disc properties was not considered in this work. 

Most of the data can be reproduced with a narrow belt of 
planetesimals centred around 40 au with sPong inward Pans¬ 
port of dust by stellar win d s, acc ording to the birth-ring sce¬ 
nario of ISPubbe & Chiangl (l2006h . Eor modelling the ALMA 
1.3 mm profile, a significant scaling factor <2 was necessary. 
This mainly hints at the inability of our simulated dust disPi- 
butions to fully reproduce the observed ALMA flux density. 
Additionally, the synthetic dust material used may underesti¬ 
mate the dust emissivity at mm wavelengths. We find a clear 
preference for a belt with a low dynamical excitation where 
the planetesimals orbits have a maximum eccentricity of 0.03. 
Our modelling allows us to make inferences on the stellar mass 
loss rate M*, which is a measure for the stellar wind strength. 
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We can exclude extreme values of hundreds of the solar mass 
loss rate, Mq, because these would cause strong deviations 
from the observed V — H and the degree of linear polarisa¬ 
tion. The best model was achieved with M* = 50 M© that cor- 
responds to the estimated wind strength during quiescent phases 
(lAugereau & Beusll2006l) . 

Due to strong dust transport from the planetesimal belt to¬ 
wards the star induced by the stellar wind drag, the inner disc 
region is filled with small scattering grains, contrary to what wa s 
derived previously (iGraham et al.lSoO^ iFitzgerald et al.ll2007ll . 
Significant deviations between the modelled and observed V-H 
and the degree of linear polarisation are visible. The problem 
may be mitigated by irregularly-shaped grains that have scatter¬ 
ing properties different from the Mie spheres used in this study. 
We find a preference for particles that are weaker polarisers than 
Mie spheres and have polarisation maxima at scattering angles 
different from 90°. This again demonstrates the necessity to in¬ 
corporate more realistic light scattering physics in simulations 
of resolved debris discs. 

Our models support the presence of ice-containing parti¬ 
cles of moderate porosity. As found by a comparison with other 
works, porous particles seem to be ubiquitous in discs around 
stars of the (3 Pic moving group. 

The radial width of the plantesimal belt cannot be con¬ 
strained tightly. Belts with radial extents of 5 and 17 au are con¬ 
sistent with the observations. Furthermore, we have addressed 
an alternative scenario with a very broad planetesimal belt, ex¬ 
tended over 1-45 au. The belt was assumed to be populated 
with particles up to 100 m in size and to have an MMSN-like 
surface density with initial radial index of -1.5. The inside-out 
evolution of this belt resulted in a rising surface density with dis¬ 
tance from the star and a 1.3 mm profile as observed by ALMA. 
The scenario explains the formation of an outer ring-like plan¬ 
etesimal distribution and would be preferred if the presence of 
planets that cleared up the inner disc region are not confirmed in 
the future. 

We have shown tha t the unresolved centra l emission at 
1.3 mm, first reported bv iMacGregor et ^ (l2013h . cannot stem 
from an inner dust ring alone. An inner ring located at < 3 au, 
which is compatible with being unresolved in the ALMA im¬ 
age, emits at most 60% of the observed central emission. We 
found a fractional luminosity ratio between inner and outer com- 
pone nt of 0.4 - a typical va lue for a two-tem perature disc sys¬ 
tem (Kennedy & Wvatdl201^ . As suggested in lMacGregor et af] 
(1201 3ll . future ALMA observations at higher resolution should 
be able to test the inner ring hypothesis. If confirmed, this might 
argue against the very extended belt scenario, based on the as¬ 
sumption that the inner planetesimals have not survived the 
inside-out collisional erosion of the disc. However, this assump¬ 
tion remains controversial because some models would explain 
the presence of an inner dust ring even after the collisional deple¬ 
tion of an extended outer ring. For example, the innermost parti¬ 
cles of the extended planetesimal belt could be less icy, having a 
higher mechanical strength, and thus, would resist the collisional 
grinding more strongly. This increases the chance to detect left¬ 
over particles from the depleted, originally extended, planetesi¬ 
mal disc. Furthermore, cometary activity may play a crucial role 
to create hot dust in the inner disc zone. Altogether, because the 
central emission cannot be explained by dust emission alone, it 
fully or at least partly derives from the stellar chromosphere. 

Our work provides a good indication of the probable generic 
architecture of the system, characterised by the spatial and size 
distributions of planetesimals and their collisionally-produced 
dust. All the models presented here are not unique as a full 


exploration of the parameter space by collisional modelling is 
impossible. However, we verified consistency of the collisional 
modelling results with a multidimensional power-law fitting ap¬ 
proach for the radial and size distributions of the dust. This re¬ 
duces the chance that our models show rather “exotic” solutions. 
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